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Disconnection of open magnetic flux by reconnection is required to balance the injection of open flux 
by CMEs and other eruptive events. Making use of recent advances in heliospheric background subtraction, 
CN '. we have imaged many abrupt disconnection events. These events produce dense plasma clouds whose 

K»" ! distinctive shape can now be traced from the corona across the inner solar system via heliospheric imaging. 

^ ' The morphology of each initial event is characteristic of magnetic reconnection across a current sheet, and 

the newly-disconnected flux takes the form of a "U"-shaped loop that moves outward, accreting coronal and 
solar wind material. 



o 
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We analyzed one such event on 2008 December 18 as it formed and accelerated at 20 m s^^ to 320 
km s^^, expanding self-similarly until it exited our field of view 1.2 AU from the Sun. From acceleration 
and photometric mass estimates we derive the coronal magnetic field strength to be 8/xr, 6Rq above the 
p ^1 photosphere, and the entrained flux to be 1.6 x 10^^ Wb (1.6 x 10^'^ Mx). We model the feature's propagation 

by balancing inferred magnetic tension force against accretion drag. This model is consistent with the 
^ ■ feature's behavior and accepted solar wind parameters. 

^ • By counting events over a 36 day window, we estimate a global event rate of Id^^ and a global solar 

minimum unsigned flux disconnection rate of 6 x lO^^Wby^^ (6 x lO^^Mxy^^)by this mechanism. That rate 
corresponds to ~ —O.lnTy^^ change in the radial heliospheric field at 1 AU, indicating that the mechanism 
is important to the heliospheric flux balance. 



1. Introduction 

The solar dynamo generates magnetic flux inside the Sun, whch is transported outward and emerges 
through the Sun's surface into the corona. Magnetic loops build up "closed" magnetic flux (connected to 
the Sun at both ends) in the corona. Some of these closed loops subsequently "open" into interplanetary 
5-J ' space - that is, they are connected to the Sun at only one end with the other extending to great distances 

in the heliosphere or beyond. Owing to the very high electrical conductivity, open magnetic flux is frozen 
into the solar wind and carried out with it. The magnetized solar wind expands continuously outward from 
the Sun in all directions, filling and inflating our heliosphere and protecting the inner solar system from the 
vast majority of galactic cosmic rays. The balance between the opening and closing of magnetic flux from 
the Sun is thus critical and fundamental both to the solar wind and to the radiation environment of our solar 
system. 

Magnetic flux opens when coronal mass ejections (CMEs) erupt through the corona, carrying pre- 
viously closed magnetic loops beyond the critical point where the solar wind exceeds the Alfven speed 
(typically <20Rq) and can no longer return to the Sun. CMEs were first studied in OS 0-7 and Skylab ob- 
servations of the corona (e.g Touseylll973 : Gosling et al iHundhausenll 1 993h . and since then continued 
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work has provided an increasingly detailed picture of these transient magnetic structures both during their 
formation and ejection, and as they continue to evolve and interact with the solar wind. 

Long lasting, radial "legs" are often observed along the flanks of a CME and persisting behind it. 
These legs are generally interpreted as evidence for at least some continued magnetic connection of CMEs 
back to the Sun and hence the opening of new magnetic flux with CME ejections. That picture is further 
supported by observatation, in situ, of beamed suprathermal halo electrons streaming in both directions 
along the lo cal interplanetary ma gnetic field (IMF) during the passage of an interplanetary CME (ICME) 
cloud (e.g. lGoslinglll990l . Il993l & references therein), which are commonly interpreted as signatures of 
direct connection of the ICME magnetic field to the solar corona in both directions, and hence of newly 
opening magnetic flux. However, less is known about the equally necessary process of disconnection that 
must be present to remove newly opened flux and prevent the IMF from growing without limit. 

Because of the continual opening of magnetic flux through CMEs, McComas and coworkers in the 
early 1990s pursued a series of studies to determine how magnetic flux could be closed back off and avoid a 
so-called r nagnetic flux "ca tastrophe" of ever increasing magnetic field strength in the interplanetary mag- 
netic field (|McComasll 1 995I & references therein). The amount of open magnetic flux in interplanetary space 
can be approximated with the "total flux integral" which removes t he effects of variations i n the solar wind 
speed in determining th e amount of m agnetic flux crossing 1 AU (|McComas et al.lll992al '). Using this in- 
tegral, McComas et al. Jl992al : ll995h showed that if all counterstreaming electron events represent simply 
connected opening magnetic loops, then for solar maximum CME rates, the amount of flux crossing 1 AU 
would double over only ~9 months. For flux rope CMEs, significantly more magnetic flux may be observed 
in the loops crossing 1 AU than what remains attached to the Sun along the CMEs' legs; however, it must 
be stressed that if CMEs retain any solar attachment whatsoever, the flux catastrophe will ultimately occur 
in the absence of some other process to close off previously open fields. 

Of course a magnetic flux catastrophe is not observed in the solar wind and, in fact, the overall mag- 
nitude of the IMF and amo unt of open flux se ems to vary over the solar cycle. For cycle 21, the average 
ma gnitude varied by ~50% dSlavin et al.lll986r) while the total flux integral varied by ~60% (McComas et 
al. 



1992a 



19920) . with maxima shortly after solar maximum and minima shortly after solar minimum 



McComaslll99 4l Since these studies, the solar wind has gone through a prolonged (multi-cycle) reduction 
in both solar wind pow er (the dynamic pressure of th e solar wind that ultim ately inflates the heliosphere) 
( McComas et al. 2008 ) and magnetic field magnitude ( Smith & Baloghll2008l) . The lack of a flux catastro- 
phe, solar cycle variation and now long-term reduction in the open magnetic flux from the Sun all show that 
there must be some process for closing off previously open field regions and returning magnetic flux to the 
Sun. 

Magnetic reconnection plays an important role in regulating the topology of solar magnetic flux, how- 
ever, once the top of a loop passes the critical point, its magnetic flux remains open until some other process 
occurs to close it off below the critical point. That is, reconnection above the critical point can only rearrange 
the topology of open magnetic flux in the heliosphere - only reconnection between two oppositely directed 
(inward and outward field) regions of open flux close to the Sun can close off previously open magnetic 
flux. The most obvious method of reducing the amount of magnetic flux open to inte rplanet ary space is via 
reconnection between oppositely directed, previously open field lines (McComas et al. Il989h . which creates 
closed field loops that can return to the Sun and the release of disconnected U-shaped field structures into 
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inteiplanetary space. An example of such a coronal disconnection event was shown by lMcComas et al. II199I 



using SMM coronagraph images from 1 June 1989. An even older exam ple of a likely coronal disconnection 
event can be found as far back as the 16 April 1893 solar eclipse (e.g.. lCliveiill989|) . where sketches (data 
in 1893) made in time or dered s equence from Chile, Brazil, and Senegal, indicate the outward motion of a 
large U-shaped structure (IMcCo mas 1994). 

For the opening and closing of the solar magnetic flux to maintain som e sort of equ ilibrium, there must 
be some type of feedback between these two processes. McComas et al. (119891 : Il99lh suggested that this 
feedback occurs through transverse magnetic pressure in the corona, where the expansion of newly opened 
field regions must enhance transverse pressure and compress already open flux elsewhere around the Sun. 
When enough pressure builds up, reconnection between oppositely direct open flux would reduce the pres- 
sure and amount of open flux. The sequence of images from the 27 June 1988 coronal disconnection event, 
in fact showed just such a compression, indicated by the deflection of the streamers in the corona, just prior 
to and appearing to precipitate the c oronal disconnecti on event. Another line of supporting evidence was 
provided by numerical simulations dLinker et al. 1992), which indicated that increased r nagnetic pressure 



could lead to reconnection across a helmet streamer and the release of disconnected flux. Schwadron et al. 



recently reexamined the flux balance issue in light of the anomalously long solar minimum between 
cycles 23 and 24 and modeled the level of magnetic flux in the inner heliosphere as a balance of that flux 
injected by CMEs, lost through disconnection, and closed flux lost through interchange reconnection near 
the Sun. 

iMcCom as et al. ( 1992c )conducted a statistical study of three months of SMM coronagraph observa- 
tions (.HundhausenI (Il993h ) to assess the frequency of coronal disconnection events. These authors found 
that while the initial survey dSt. Cyr & Burkepildl 19901) found no obvious disconnections, six of the 53 tran- 
sient events during this interval (11%) showed some evidence of disconnection in more than one frame and 
13 (23%) showed a single frame with an outward "U" or "V" structure. Given the imaging and analysis 
technology of the day, McComas et al. (1992c) concluded that magnetic disconnection events on previ- 
ously open field lines may be far more common than previously appreciated. With today's imaging and 
exceptional analysis capabilities, the question of coronal disconnection events should finally be resolvable. 



For this study we used image sequences, collected by the SECCHI (IHoward et al.ll2008 r) instrument 
suite on board NASA's STEREO-A spacecraft, of Thomson-scattered sunlight from free electrons in the in- 
terplanetary plasma. The observations span from the deep solar corona to beyond 1 AU at elongation angles 
of up to 70° from the s olar disk; this continu ous observation is enabled by recently developed background 
subtraction techniques ( DeForest et al.ll2011 ) operating on the STEREO data. The signature U-shaped loops 
of disconnected plasma are far clearer in the processed heliospheric data far from the Sun than in the coro- 
nagraph data close to the Sun, and we detect 12 characteristic departing "V" or "U" events in 36 days - 
far more than the expected number based on scaling the results of McComas et al. (1992c). For this initial 
report, we focus on quantitative analysis of a single event. In Sections 2.1-2.7, we describe the observations 
and calculate the geometry, the mass evolution, and (by assuming the U-loop is accelerated by the tension 
force) the coronal magnetic field sand entrained flux in the disconnecting structure. As a plausibility check, 
we explore the tension force scenario and its consequences for the long-term evolution of the feature, and 
find that the scenario is consistent with accepted values for the solar wind density and speed. In Section 3, 
we discuss broader consequences of the observation, including estimating the disconnection rate based on 
the number of similar events in our data set, and discuss implications for the global magnetic flux balance. 
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2. Observations 



The SECCHI s uite on STEREO was intended to be used as a single integrated imaging instrument (e.g. 
Howard et al.ll2008r) . It consists of an EUV imager (EUVI) observing the disk of the Sun, and four visible 
light imagers (COR-1, COR-2, HI-1, and HI-2) with progressively wider overlapping fields of view, to cover 
the entire range of angles between the solar disk and the Earth. The visible light imagers view sunlight that 
has been Thomson scattered off of free electrons in the coro na and interplanetary space ; the theory of Thom- 



son scattering observations has been recently reviewed by iHoward & Tappin 



2009al ). We set out to 



view 



coronal and heliospheric events in the weeks around 2008 December, using newly develop ed background 



subtra c tion techniques to observe s olar wind features in the HI-1 and HI-2 fields of view (IDeForest et al. 



(120 lib : Howard & DeForesli (1201 ih ). In the initial 36 day data set we prepared, we observed 12 disconnec- 
tion events identified by a clear "V" or "U" shaped bright structure propagating outward in the heliosphere. 
We chose a particuarly clearly presented one, which was easily traceable to its origin in the low corona on 
2008-Dec-12 at 04:00, for further detailed study. 



2.1. Image Preparation 

Data preparation followed standard and published techniques. For COR-1 and COR-2, we downloaded 
Level 1 (photometrically calibrated) data, and further processed them by fixed background subtraction: we 
acquired images for an 11 -day period, and found the 10 percentile value of each pixel across the entire 
11 -day dataset. This image was median filtered over a 5x5 pixel window to generate a background image 
that included the F corona, any instrumental stray light, and the smooth, steady portion of the K corona. 
Subtracting this background from each image yielded familiar coronal images of excess feature brightness 
compared to the smooth, steady background. We performed one additional step: motion filtration to suppress 
stationary image components. This step matches the motion filtration step used for HI-1 and HI-2 (below), 
and suppresses the stationary streamer belt while not greatly affecting the moving features under study. 

The heliospheric imagers required further processing to remove the starfield, which is quite bright 
compared to the faint Thomson scatt ering signal far from t he Sun in the image plane. We processed the 



STEREO-A HI-2 data as described by iDeForest et al.l (|201 11 ). The HI-1 data used a similar process a dapted 



to the higher background gradients in that field of view and described by lHoward & DeForesti (|201ll ). 



All the imagers yielded calibrated brightness data in physical units of the mean solar surface brightness 
( 50 = 2.3 X lO^W m^'^ SR^^). Because of the wide field of view, as a subsequent processing step we 
distorted the images into azimuthal coordinates, in which one coordinate is azimuth (solar position angle) in 
the image plane and the other is either elongation angle e ("radius" on the celestial sphere) or its logarithm. 
The latter projection, if scaled properly, is conformal: it preserves the shape of features that are small 
compared to their distance from the Sun. To equalize brightness, we applied radial filters to the images for 
presentation, with either a e^'^ scaling (for coronal images) or a scaling (for heliospheric images). 

Figure[I]shows snapshots of the disconnected plasma and associated cusp, as observed by four separate 
instruments over the course of four days as it propagated outward. Shortly after 2008 Dec 18 04:00, the 
streamer belt at 160° ecliptic azimuth (20° CCW of the Sun-Earth line) pinched and separated, forming a 
"U" loop that retracted outward, with a trailing cusp, over the course of the following three days. The feature 
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remained visible in Thomson scattered light because of plasma scooped up during the early acceleration 
period in the lower corona: this plasma remained denser than the surrounding medium, yielding a bright 
feature throughout the data set. The disconnected plasma completely missed the ecliptic plane and was 
therefore not observed in situ by any of the near-Earth or STEREO probes. 



2.2. Observing Geometry and 3-D structure 



The observing geometry for the 2008 Dec 18 event is shown in Figure |2l from an overhead (northward 
out-of-ecliptic) point of view co-rotating with the STEREO-A orbit. The event departure angle was mea- 
sured both using direct triangulation between the coronagraph s in STEREO-A and STEREO-B. We used 
the triangulation method described by iHoward & TappinI (|2008r) . Although the disconnection event is small 
compared to most CMEs, subtending just a few degrees in latitude, it is still large enough to cast doubt 
on the si mple triangula tion results, so we also used "TH model" semi empirical transient event reconstruc- 
tion tool (^Ta ppin & Howard.200 9) to extract the departure angle. TH was developed to reconstruct CME 
leading edge ("sheath") overall envelope and propagation speed, but is also applicable to smaller transient 
events such as this one . Details, applications, an d limitations of the TH model are further described by 
Tappin & Howardll2009l and by Howard & Tappin (l2009bl : l201Clh . Departure longitude was measured to be 
— 10°+5° in heliographic coordinates, with an estimated event width of under 5°. 

We took the disconnected feature's trajectory to have constant radial motion (the "Fixed-O aproxima- 
tion") in the solar inertial frame - this leads to the slightly curved aspect to the trajectory in the co-rotating 
heliographic ecliptic frame, which maintains the prime meridian at the Earth-Sun line. Figure [2] shows an 
out-of-ecliptic projected view of the observing geometry, including construction angles and distances used 
in Section l24l for trajectory calculations. 



2.3. Feature Evolution 



To analyze the feature's evolution across a two-order-of-magnitude shift in scale over its observed life- 
time, we transformed the processed STEREO-A source images into local heliographic radial coordinates - 
i.e. zero azimuth is due solar West from the viewpoint of STERE O-A, with azimuthal c oordinate increas- 
ing clockwise around the image plane; this follows early work by iDeForest et all (120011) in imaging polar 
plumes. Distances from Sun center are recorded as elongation angle e from the center of the Sun, as a 
reminder of the angular nature of the wide-field observations. To avoid aliasin g in the re s ampli ng process, 
we resampled the images using the optimized resampling package described by lDeForesti (120041) . Figures [3] 
and |4] show the liftoff and propagation of the feature across 65 degrees of elongation from its origin in the 
solar streamer belt. Both figures have a radial gain filter applied to equalize the feature's brightness, which 
varies by over seven orders of magnitude: from 1.5 x 1O^^S0 in the low streamer belt at 2008 Dec 18 04:30 
to 5.7 X 1O"^''S0 six days later, at £ = 65°. 

The bright feature t akes the classic w ishbone shape of reconnecting field lines emerging from a current 
sheet (e.g. Chapter 4 of IPriest et al.ll200Clh . The aspect ratio of the wishbone may be estimated by dividing 
the vertical height from cusp to the top of the visible horns, by the width between the horns. This aspect 
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ratio varies from ~10:1 when the horns are first clearly resolved near 2008 Dec 18 08:00, to approximately 
2:1 some four hours later and 1:1 by 2008 Dec 19 04:00 - one full day after the first pinch is observed in 
the streamer belt. After 2008 Dec 19, the feature expands approximately self-similarly as it propagates, sub- 
tending approximately 16° of azimuth and not changing its aspect ratio throughout the rest of its trajectory. 

Note that aspect ratio is not preserved by the linear azimuthal mapping used in Figure [3l which was 
selected to show the early acceleration clearly; aspect ratio is preserved by the (conformal) logarithmic 
mapping used in Figure |4l which shows nearly self-similar expansion in the image plane despite perspective 
effects that come into play above about e = 30° 

The scaling of brightness is reassuring because, in a uniformly propagating wind with no acceleration, 
density must decrease as r^^ and feature column density must thus decrease as r^', while illumination 
decreases as r^^, so feature brightness is expected to decrease as r^^. The fact that brightness levels do 
not change much across Figured which is scaled by suggests that the disconnected flux and material 
entrained in it are indeed propagating approximately uniformly. The fact that they do change slightly, with 
brighter images to the right, indicates that the feature is gaining intrinsic brightnsss by accumulating material 
as it propagates. 

The horizontal positions and error bars in Figures [3] and |4] are the results of manual feature location of 
the cusp, with a point-and-click interface. The white error bars are based on the sharpness of the feature. 
In the excess brightness plot, the feature is easy to see but blurs near the top of Figure [3] due to the higher 
levels of both photon noise and motion blur as the feature accelerates to the top of the coronagraph field of 
view. The running difference plot highlights fine scale feature and helps identify the cusp location near the 
top of the COR-2 field of view. 



2.4. Acceleration profile 

Converting angular observed coordinates to examine the inertial behavior of the plasma requires trian- 
gulation using the Law of Sines. Using the "fixed <I>" approximation (assuming the feature's cusp is small 
and that it propagates in a radial line from the Sun), the feature's radius from the Sun is easily calculated: 

sin{e') 

rev = rA . , (1) 
sin (e) 

where the variables take the meanings in Figure |2l £ is the solar elongation of the feature as seen from 
STEREO-A, e' is the solar elongation of STEREO-A as seen from the feature), and the STEREO-A solar 
distance is found by spacecraft tracking and is supplied by the mission. Although no camera was present 
at the event itself, e' is calculated by noting that = 180° — £ — (L — Lev), as the feature, STEREO-A, and 
the Sun form a triangle. Figure [5] shows the results of the tracking from Figures [3] and IH propagated through 
Equation [T] 

As expected, the event rapidly accelerates during the early phase, reaching a peak acceleration of 
20ms^^ as the aspect ratio changes in the initial hours. The acceleration peaks 4-5 hours after the initial 
pinch in the streamer belt, or 2-3 hours after the first observation of a well-formed cusp. The feature reaches 
its final speed of ~ 320± I5kms^^ within just 8 hours of the initial pinch at 04:00 and within 6 hours of 
the first observation of the well formed cusp at 06:00, and undergoes no further significant acceleration nor 
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deceleration during its obsered passage to beyond 1 AU over the next five days. 

2.5. Mass profile 

We extracted photometric densities using the feature brightness in suitabie frames. The feature bright- 
ness is determined from the density via the Thomson scattering equation (see, e.g., Howard & Tappin 2009a 
for a clear exposition). Compact features can be treated as nearly point sources, and the line-of-sight integral 
for the optically thin medium reduces to: 

B = BQaQ{r)ae{l+cos^x)p^tavd (2) 

where B is the measured feature brightness (in units of emissivity: Wm^^SR^^), Bq is (still) the solar 
surface brightness; ^Q{r) is the solid angle subtended by the Sun at the point of scatter, well approximated 
by Tlr^r^^ everywhere above about A r^; is the differential Thomson scattering cross section, given by 
half of the square of the classical electron radius rl/2 = 4.0 x 10^^°m^; % is the scattering angle (equal to 
e' in Figure O; p is the mass density; }Xav is the average mass per electron in the coronal plasma; and d is 
the depth of the feature. 

muav may be calculated from the spectroscopically measured 5% He/H number ratio in the corona 
(Laming & Feldman 2000) and the assumption that the helium is fully ionized (yielding two electrons per 
ion). This yields mUav = l-lnip = 1.84 x 10^^^ kg. 

Solving for the line-of-sight integrated mass surface density pd gives: 

pd = ^av^Q.-\r)o-' {l+coshy' (3) 

Bq 

and therefore 

m,, = {pd)wh = ^^,^^-\r)a-' {l+cos^e'y'dl^""'''h^ (4) 

Bq 

where w and h are the dimensions shown in Figure |7J D.ev is the solid angle subtended by the feature in the 
images; and 5 is the calculated spacecraft-feature distance, calculated by the law of sines as for rev- 
To extract the mass profile from the data, we generated an image sequence containing the feature, and 
marked the locus of the feature visually using a pixel paint program. Using the generated masks, we summed 
masked pixels in the feature for each photometric image, thereby integrating the feature brightness over the 
solid angle represented by the corresponding pixels, to obtain an intensity and an average brightness within 
the feature. To account for errors in visual masking, we assigned error bars based on one-pixel dilation and 
one-pixel contraction of the masked locus. We omitted frames with excessive noise, encroachment of an 
image boundary, or a star or cosmic ray in or near the feature. The results of the calculation are given in 
Figure [6l which shows steady accretion of material through most of the journey through the heliosphere. 

Because our photometric analysis is based on subtraction of a calculated background derived from the 
data set itself, we measure only excess feature brightness (not absolute brightness) from Thomson scatter- 
ing; thus our brightness measurements and mass estimates are biased low, because we cannot measure the 
absolute density of the background. The initial derived mass of 20-25 Tg translates to an electron number 
density of 2 x 10^ cm^^ in the lower corona, which is comparable to the density in bright coronal features - 
so the total mass may be up to a factor of order two higher. 
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The final "feature excess" mass is 8 ± 2 x 10 kg, and the final subtended solid angle is 0.028 SR, for 



a presented cross-section of 4.3 ±0.1 x 10 m . Taking the depth to be the square root of the observed 
cross section yields an estimated volume at 1 AU of 8.9±0.3 x 10^'m^, for a total estimated excess electron 
density of 5 ± lAan^^ at 1 AU, which is in good agreement with slow solar wind densities (3 — lOcm^^ 
when scaled to 1 AU ) that were observed by Ulysses in situ in the same heliographic latitude range (e.g. 



McComas et al.ll2000l '). Approximately 2/3 of this excess density appears to have been accumulated enroute 



from the surroundign solar wind; this is further described in Section 12771 



2.6. Entrained magnetic flux 

From the mass of the feature, and its acceleration, it is possible to extract the entrained magnetic field 
by measuring the rate of change of momentum and inferring a magnetic tension force via / = ma. The 
system is sketched in Figure |7] The magnetic tension force is conserved along the open field lines, so we 
can calculate it at any convenient cut plane including the one shown. Tension force is frequently referred to 
as a "curvature force" and calculated locally; here we integrate around the "U", and notice that the integrated 
force is just the unbalanced tension on the field lines contained in the "U" shape. It is therefore given by 

/ = m,„..„ = /, = — = ^ (5) 
where, here, B is the magnetic field strength (not brightness, as before). Solving for 0, 



<I> = \/ lHodwrrievaev (6) 

taking nigy to be 25Tg (2.5 x lO^^kg) during the peak of the acceleration, and taking w = d = 0.2Rq (based 
on the measured width of the feature's fork during maximum acceleration, at 6Rs from the surface (7/?^ 
from Sun center) gives = 1.6 x lO^^Wb (1.6 x lO^'^Mx ), corresponding to an average field strength of 
8/1 r ( 0.08 Gauss) at that altitude, or an equivalent r^-scaled field of 400/1 T (4 Gauss) at the surface; this 
is comparable to accepted values of the open flux density at the solar surface at solar minimum. Because of 
the way m was calculated (section 2.5 above) this figure is probably low by a factor of order \/2. 



2.7. Accretion and force balance 

As the disconnectioned structure travels outward, it accretes new material. This effect is dramatic: as 
seen in Figure[6l the mass increases by a factor of 3 from the corona to 1 AU. We conjecture that the material 
is accreted by "snowplow" effects from the plasma ahead of the disconnected cusp as it propagates. For the 
observed mass growth in the feature, new material must be compressed to become visible in our Thomson 
scattering images, and the most plausible way for it to be compressed is via ram effects. This scenario also 
neatly explains the constant speed of the feature, by balancing the continued tension force from the cusp 
with accretion momentum transfer. Here we explore the concept of force balance between accretion and the 
tension force, to identify whether some other model is required in addition to this simple one. 

Extending Newton's law to include momentum transfer by accretion, and neglecting all but the tension 
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force, 



<J)2 



dm, 



■Av, 



2/ioA<i> 



(V) 



where the LHS is just the tension force from Equation |5l with the modification that the cross section of 
the exiting field lines is written A^; = after the initial acceleration; and the second term represents 
momentum transfer into accreted material, with Av being the difference between the feature speed and sur- 
rounding wind speed. The feature is thus in equilibrium between accretion drag and continued acceleration 
by the tension force. This accretion drag is important to the observed increase in feature mass, because ram 
pressure against the surrounding wind material is what compresses incoming material and renders it visible 
in the data against the subtracted background. 

Applying conservation of mass, we can relate Av and the average density of the background solar wind 
through which the feature is propagating: 



where Ag,, is the geometrical area presented by the feature to the slow wind ahead of it. Solving Equations |7] 
and [8] to eliminate Av gives 



which gives the background solar wind density in terms of the accumulation rate of mass in the observed 
feature, assuming constant outflow for both the wind and the feature, and acceleration by the tension force. 
Given the conservation of mass and the approximately constant speed of the solar wind, Psw falls as ap- 
proximately r^^. Further, we observe nearly self-similar expansion throughout most of the heliospheric 
range, so Aci)/Aev is constant in that part of the trajectory - hence dmev/dt must also fall as r^^ during the 
approximately constant speed portion of the feature's lifetime. Using this functional form, we can extract an 
analytic expresson for the feature mass versus radius. We introduce the r^^ dependence by switching from 
the linear regression used in Figure [6l to a semi-log regression that assumes dmev/d{loge{r) to be constant. 
Figure [8] shows such a regression, with the result that diriev/ <ir = 34 it 3 x 10^ {Rq / r) kgR^ ' . Including the 
measured outflow speed of 315 it \5kms^^ , we find that dmev/dt = {\AU jr) (7.1 it 1 x 10'*%5'~^). 

Including all of these values into Equation |9j together with the average particle mass from Section [231 
yields a background wind numeric density of ?i„, (lA?7) of 30ib6cm^^ (A^/Agv). From the morphology 
of the feature in Figure |4j we conservatively estimate A^/A^v < 0.25, i.e. the forward cross section of the 
"horns" of the vee appears to be well under 1/4 of the cross section of the vee itself. This value yields a 
derived background solar wind density of ?is„ < 8cm^^ at 1 AU to maintain the force balance in Equation |7] 
That figure is again in line with the wind measurements from Ulysses at 15° heliographic latitude (McComas 
et al. 2000), adding to the plausibility of the accretion force balance picture. The corresponding mass density 
limit is p^v, < 8 X 10^^'^ kg m^^ at 1 AU 

As a sanity check, we can use this p^v limit and Equation [8] to find that Av must then be a few tens of 
kms^^, i.e. the background wind speed must be close to 300kms^^. 

We conclude that the picture of force balance between snowplow accretion and the tension force is 
at least broadly consistent with the observed feature, though further study of more events (preferably with 
corresponding in situ measurements of the feature itself) is necessary. 



Psw 



dmev/dt 



(8) 



A«.,,Av 




(9) 
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3. Discussion 

Using data from STEREO/SECCHI, we have identified and measured the characteristics of a single flux 
disconnection event and associated cusp feature, similar to that discovered by McComas et al. (1992), from 
initial detection in the lower corona to distances beyond 1 AU. The cusp feature is formed in the classic X- 
point geometry and rapidly accelerates under the tension force to approximately 320 km s^^, which it reaches 
in under 4 hours at an altitude of approximately IORq. Thereafter the feature continues to accumulate mass 
but maintains approximately constant speed until it is lost to sight 1.2 AU from the Sun. 

Based on photometry, we are able to estimate the onset mass of the event as 25 Tg and the entrained 
flux as l60GWb, corresponding to a coronal field strength of 0.08 G and an -normalized surface open 
field of 4 G over the projected surface footprint of the feature. These estimates are Ukely low by a factor 
of order \/2, because they make use of feature excess brightness rather than absolute Thomson-scattered 
brightness in the coronagraph images; using polarized-brightness imagery could improve the measurement 
by separating the non-transient component of the Thomson scattering signal from the unwanted F coronal 
background. 

Because our measurements are all based on morphology and photometry, we have performed several 
consistency checks to build confidence in the calculated parameters of the feature as it propagates. In 
particular, a model of simple force balance between the tension force and mass accretion is consistent with 
both the inferred magnetic field and accepted values for background slow solar wind density and speed. 

Simple accretion models such as we developed here demonstrate clearly why ejected features such as 
U-loops or CMEs seem frequently to propagate at near constant speed: under continuous weak driving, 
an equilibrium forms rapidly between the driving force and momentum transfer by mass accretion. The 
equilibrium outflow speed is the sum of a large, fixed (or at least driver-independent) speed - that of the 
surrounding wind - with a smaller offset speed that drives mass accretion. Thus the feature speed is quite 
insensitive to the driver. In our case, doubling the tension force would only increase the outflow rate by 
~ 10%. 

The event under study is well presented, but is not unusual at all; such events are easy to identify in 
heUospheric image sequences, because of their distinctive "U" and cusp shape; they are readily traced back 
to the corona. This technique represents a new, very effective way of finding these discoimection events, 
which are small and hard to identify in the coronagraph sequences alone, but are strongly and easily visible 
in the processed heUospheric images. 

In an initial reduced data set of 36 days near the deepest part of the recent extended solar minimum 
(2008 Dec - 2009 Jan), we identified 12 such events; all of them were identified by tracking "V" or "U" 

shapes back from the heliospheric images to the corona. Assuming the present feature to be typical, and 
considering that the single viewpoint affords clear coverage of about 1/4 of the circumference of the Sun, we 
estimate the global disconnection feature rate at that time to be over I event d^^, and the flux disconnection 
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rate to thus be at least of order SOTWby^ . Expanded to a 1 AU sphere, this amounts to a rate of change 
of the open field of order O.lnTy^^, which is a significant fraction of the observed cycle-dependent rate of 
change of the open heliospheric field (e.g. Schwadron, Connick, & Smith 2010). These figures are based 
on a single calculated flux and an event rate obtained by initial visual inspection of a single 36-day data set, 
and hence are merely rough estimates - but they indicate that flux disconnections of this type are important 
to the global balance of open flux. Further study, in the form of a systematic survey, is needed to determine 
whether they are the primary mechanism of flux disconnection from the Sun. 

The authors thank the STEREO instrument teams for making their data available. Our image process- 
ing made heavy use of the freeware Perl Data Language (http://pdl.perl.org). The work was enhanced by 
enlightening conversations with J. Burkepile, C. Eyles, and N. Schwadron, to whom we are indebted. This 
work was supported by NASA's SHP-GI program, under grant NNG05GK14G. 
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Fig. 1. — The disconnection event of 2008 Dec 18 in context - 8 still images showing formation and evolution 
of the U-loop: left to right, top to bottom. 



-14- 




mm imo mm io:oo 12:00 14:00 wm 




08:00 10:00 1» 1« 



Fig. 3. — Formation and early acceleration of the 2008 Dec 18 disconnection event through the STEREO- 
A COR-1 and COR-2 fields of view The traiUng edge of the event is marked, with error bars based on 
feature identification. TOP: direct excess-brightness images show feature formation and overall structure. 
BOTTOM: running-difference images show detail. These stack plots include a small image of the feature at 
each sampled time to show evolution. Intensities are scaled with e^-^ to equalize brightness vs. height. The 
individual images have been resampled into linear azimuthal (radial) coordinates, and the horizontal range 
is 160° — 174° of azimuth. Note that this projection does not preserve aspect ratio: despite appearances, the 
event widens as it rises. 
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Fig. 4. — Propagation and evolution of the 2008 Dec 18 disconnection event through the STEREO- A HI-1 
and HI-2 fields of view. The trailing edge of the event is marked, with error bars based on feature identifi- 
cation. These stack plots include a small image of the feature at each sampled time to show evolution. The 
individual images have been resampled into logarithmic azimuthal (radial) coordinates, and the horizontal 
range is 162° — 178°of azimuth. This projection is conformal, so the shape of the feature is preserved in 
each image. Intensities are scaled with to equalize brightness vs. height. Note self-similar expansion: the 
angular width and shape of the feature are preserved. 
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Fig. 5. — Inferred position, speed, and acceleration of the disconnected plasma from the 2008-Dec-18 event, 
during onset (LEFT) and over the full observation period (RIGHT). Error bars are derived by propagating a 
priori location error and geometric error in the longitude of the event. The shaded region indicates the full 
time range of the left-side plots. 
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Fig. 6. — Photometrically determined excess mass profile of the retracting disconnected feature of 2008 Dec 
18. Error bars are based on identification of the feature boundary in the images. The trendline is extracted 
from regression of the HI-1 and HI-2 data. The mass shown is excess mass in the feature compared to the 
background solar wind (see text). 
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Fig. 7. — Cartoon of the initial acceleration process of a disconnection event. Tension force along newly 
released field lines is balanced by mass entrained on the field fines. By measuring the acceleration and mass 
we infer the amount of magnetic flux that was disconnected. 
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Fig. 8. — Semilog regression fit of mev{r) 



